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functions 



Markus Spanner 3 , Simon K. Schnyder b , Felix Hofling c , Thomas Voigtmann d , and Thomas Franosch 

Received: November 20, 2012 



o 

(N 

> 
O 

On 



o 

S 

i 



I 

-a 
c 

o 
o 



> 
o 
m 
in 



Heterogeneous media constitute random disordered environments where transport is drastically hindered. Employing extensive molecular 
dynamics simulations, we investigate the spatio-temporal dynamics of tracer particles in the Lorentz model in the vicinity of the localization 
transition. There transport becomes anomalous and non-gaussian due to the presence of self-similar spatial structures, and dynamic scaling 
behavior is anticipated. The interplay of different time and length scales is revealed by the intermediate scattering functions, which are 
sensitive both to the underlying spatial fractal as well as the anomalous transport. We compare our numerical results in the transition regime 
to a mode-coupling approach, and find that certain aspects are surprisingly well predicted. 
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1 Introduction 

One of the most striking phenomena of densely packed matter is 
the rapid slowing down of dynamics upon cooling or compression 
as the glass transition is reached. Certainly, the most prominent 
aspect is the increase of viscosity by orders of magnitude " upon 
changing the temperature within a few degrees or similarly the 
density by a few percent. At the same time the diffusion of the 
particles is strongly suppressed suggesting that motion is strongly 
hindered until essentially a structural arrest occurs. More insight 
on the evolution of the glass transition is obtained by focusing at 
time- or frequency-dependent quantities that encode the dynamic 
processes and the system's response on different temporal scales^] 
as can be nicely achieved in colloidal systems.El The next step is 
to probe the sample also at different length scales to unravel the 
mechanism that is responsible for the drastic slowing down of dy- 
namics. A series of non-trivial phenomena such as the stretching 
of time-dependent correlation functions, the time-temperature su- 
perposition principle, the emergence of a fast /j-process, where a 
factorization property of the space and time dependence occurs, 
have been observed in the last 25 years. 2 3 A coherent theoretical 
framework has been provided in terms of the mode-coupling theory 
(MCT) of the glass transition developed by Gotze and collabora- 
tors.^ Although many predictions of this microscopic theory have 
been tested successfully in experiments and computer simulations, 
the limitations of the MCT approach are poorly understood. 

Recent extensions of MCT^nSJto liquids confined inside a porous 
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medium predict a rich interplay of the glass-transition singu- 
larity and the localization phenomenon in disordered structures, 
parts of which have been verified already by computer simula- 
tions EHE] T nere me p 0rous environment serves as a frozen host 
structure that strongly interacts with the dense liquid filling the void 
space. This scenario appears to be realized approximately in ion- 
conducting sodium silicate melts,! 18 " 20 l where the silicon-oxygen 
majority component displays diffusion coefficients that are by or- 
ders of magnitude smaller than the ones of alkali ions. Thus the 
alkali ions essentially meander through an arrested structure and 
move along certain preferential diffusion pathways percolating the 
matrix. The scenario of a multiple freezing transition has been 
also discussed for strongly size-disparate colloidal mixtures of soft 
spheres 21 " 24 or Yukawa particles. 25 

In the limit of very small concentration of the liquid inside the 
porous structure, the system coincides with the Lorentz model, 
where a gas of non-interacting tracer particles explores a frozen 
disordered array of obstacles. There the transport properties de- 
pend sensitively on the obstacle density and eventually long-range 
transport ceases to exist as a certain critical density is reached. Ex- 
tensive computer 

simulations^!] 

confirm that the mean-square 
displacement increases subdiffusively directly at the transition and 
scaling behavior holds in the close vicinity of the critical density. 
A self-consistent mode-coupling kinetic theoryEHS] ma de a series 
of predictions of the dynamic behavior for all obstacle densities, 
in particular, it reproduces the long-time anomalies in the veloc- 
ity autocorrelation 

functions£212U 

and yields an accurate estimate 
for the first-order correction. 35 Furthermore the theory predicts a 
transition to a localized phase and suggests spatio-temporal scal- 
ing behavior in its vicinity. The mechanism of the suppression of 
diffusion is due to repeated correlated scattering with the obstacles 
confining the tracer to the still accessible void space. Yet, in the 
immediate vicinity of the transition a new phenomenon becomes 
important. The excluded volume due to obstacles starts to overlap 
significantly and builds clusters, eventually leading to self-similar 
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structures for the void space. Hence a geometric percolation tran- 
sition gives rise to a divergent static length scale which ultimately 
changes the critical behavior. 

The purpose of this paper is to explore at the transition and to 
characterize transport properties in terms of the intermediate scat- 
tering function for particles located on the infinite cluster only. 
There processes on different length scale can be monitored by 
adjusting the wavenumber. We elaborate a series of scaling re- 
lations at the transition both within a simplified mode-coupling 
approacrJSD and dynamic scaling hypothesis^ for the universal- 
ity class of the critical properties or random resistor networks. We 
compare both sets of predictions to extensive computer simulations 
and discuss the scaling properties. 

2 Lorentz model 

2.1 Model definition 

The Lorentz model is a simple model for transport in porous mate- 
rials, where a point-like tracer moves through an array of quenched 
spherical obstacles. Introduced by H. A. Lorentz 38 in 1905 to de- 
scribe the electron conductance in metals, it is used more generally 
for transport phenomena in disordered systems, and even serves as 
a microscopic model for anomalous transport in crowded biologi- 
cal media.1221 Obstacles are distributed independently in the sample, 
hence the statistical properties of the structure are simply charac- 
terized by the number density of obstacles n = N/V. Here V = Lr 
denotes the volume of a hypercubic box, ultimately we are inter- 
ested in the limit of large system sizes, N — > °°, V — > °°, while 
keeping n = const. 

The disordered environment is explored by a single structureless 
particle, equivalent to a gas of non-interacting particles. We con- 
sider a hard-core exclusion of radius cr between a single obstacle 
and the tracer. The interaction distance cr can be interpreted as the 
sum of tracer and obstacle radius, which include the limiting case 
of an extended tracer in point-like scatterers, as well as a point-like 
particle meandering in a parcours of possibly overlapping spheres. 
Hence the dimensionless control parameter is the reduced obstacle 
density n* = na^ , where d is the spatial dimension. The particle's 
motion is confined to the void space and one anticipates already 
that with increasing obstacle density the void space decomposes 
into a collection of smaller and smaller clusters. Above the per- 
colation transition n* > n* there is no cluster that spans the entire 
system and long-range transport ceases to exist. 

For the tracer dynamics, we employ Newton's equations of mo- 
tion with a specular scattering at the obstacles. Then the kinetic en- 
ergy of the tracer is conserved and only the direction of the velocity 
v is changed in scattering events, such that the velocity v = |v| sets 
the time scale t B = cr/v of the problem. 

2.2 Scaling theory of the critical dynamics 

At the critical obstacle density n*, the infinite cluster becomes self- 
similar 40 characterized by a fractal dimension df, which evaluates 
approximately to^ df = 2.53 in 3 dimensions. The dynamics on 
this fractal is expected to become anomalous since the tracer has 
to explore a network of ramified structures. The simplest quan- 



tity probing the transport properties is the mean-square displace- 
ment c>r 2 (f) := (A/JI/) 2 )^, where (•}«> indicates averaging differ- 
ent realizations for particles moving on the infinite cluster only and 
AR(f) = R (?) — R(0) is the displacement of the tracer after a lag 
time t. 

Since the percolating void space is scale-free one anticipates 
8r^(t) ~ ; 2 / rf »- 5 where <i w is known as the walk dimension of 
the system.^ F r the three-dimensional Lorentz model the value 
d w = 4.81 differs from the one in random resistor networks on a 
lattice, since transport is dominated by narrow channels emerg- 
ing in continuum percolation.l^SI^The mean-square displacement 
is connected to a corresponding time-dependent diffusion coeffi- 
cient Doo(f) := (l/2d)d(5r 2 (r)/df, and a velocity auto-correlation 
function Z„(t) := (l/2d)d 2 (5r 2 (r)/df 2 . In the Fourier domain, the 
second moment can be transformed to a frequency-dependent con- 
ductivity which displays dispersive transport directly at the critical 
density. 30 Few is known beyond the second moment, only recently 
spatial-temporal properties have been discussed.^ 

A complete characterization of the statistical properties is given 
by the van Hove function^ Poo(r,t) := (6{AR(t) - r)}oo, which 
constitutes the probability for the particle to have traversed a dis- 
tance r in lag time t. In particular, the mean-square displacement 
represents just the second moment of P„(r,t). 

Self-similarity suggests a scaling 

behaviorEzEH 

PooM = r rf ^( r r" , '«) (i) 

for large distances r 2> a and long times / 2> t . Thus one expects 
typical excursions of the tracer of linear extent ~ t , larger dis- 
tances become increasingly rare, suggesting ZP^x S> 1) to decay 
rapidly. 

The behavior for small rescaled distance x ~ rf~ ^ is in- 
ferred by the following consideration: 45 The return probability 
n(f,w) := / r < w d .rPa,(r,i) represents the likelihood that the parti- 
cle remains or has come back within a distance w from the starting 
point after a time t. Provided that this distance is much smaller 
than the typical excursions, w <C f'^ w , this probability is propor- 
tional to the accessible volume ~ w^ 1 of the infinite cluster. By the 
scaling law, eq. |TJ, this implies 8P^{x <C 1) ~ x df . 

The intermediate scattering function (ISF) 

F.( 9 ,0:=(«ap(i«.A«(0)>- (2) 

is the characteristic function of the displacement of the tracer 
AR(t) after a lag time t. This quantity is in principle experi- 
mentally accessible by neutron spin-echo spectroscopy or photon- 
correlation spectroscopy, where hq is identified with the momen- 
tum exchanged between the neutron/photon and the sample] 4 ^ Ob- 
serving that the ISF represents merely the Fourier transform of the 
van Hove function, Foo(q,t) = Jexp(i^- r)P^(r,t) d d r, one can de- 
rive a scaling relation for the critical behavior. Hence from eq. |TJ 
one concludes 

F„( q ,t) = ^(qt 1 ' d «) (3) 

for small wavenumbers q <C cr~ 1 and long times (>f„. Substitut- 
ing r i — y xt l / d ™ and q n» Kt~ l l d '" in the Fourier transform yields the 
scaling function as ^„ (k) = J exp (i*c • je) (x)d d x. In particular, 
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the return probability argument determines the large-wavenumber 
behavior &„{k > 1) ~ K~ d < 

In the vicinity of the critical density scaling still holds, how- 
ever, the divergent correlation length £ ~ (n* — n*)~ v has to be 
introduced.!^ The approach towards the scaling law is governed 
by universal corrections to scaling, which have been worked out 
recently.^ 

2.3 Mode-coupling theory 

A theory based on first principles has been developed in terms of 
a self-consistent mode-coupling kinetic fheoryEEH f or me entire 
range of densities. Here we focus on the critical regime to ratio- 
nalize the anticipated scaling behavior. The universal aspects are 
captured by a simplified approach 36 which considers the system 
composed of a tracer and the porous host structure as a binary mix- 
ture, where one component is frozen permanently. The equations 
of motion then readily follow by specializing the mode-coupling 
theory of the glass transition to the case considered here. 

First, the intermediate scattering function fulfills the exact equa- 
tion of motion 



P{q,t)+v s F(q,t)+Cl{q) 2 F{q,t)- 



- f'M(q, 
JO 



■t')F(q t t')dt' = 
(4) 



as can be derived within the Zwanzig-Mori projection operator 
formalism.^ For definiteness we specify the discussion in the re- 
mainder of the paper to d = 3. The frequency Q.{q) 2 = q 2 v 2 /3 
corresponds to the ballistic motion of the particle. The interac- 
tion with the disordered matrix leads to memory effects and fric- 
tion. The collision rate v s = nncr 2 v induces a Markovian damp- 
ing and accounts for uncorrelated scattering with the host matrix, 
while the non-trivial correlations are hidden in the memory ker- 
nel M(q,t). Within the mode-coupling approximation, the mem- 
ory kernel is split into a regular smoothly changing background 
and a mode-coupling contribution that entails a feed-back mech- 
anism. For simplicity we keep only the latter term, bearing in 
mind that the Markovian damping should be adjusted to account 
for the smoothly changing background. The mode-coupling kernel 
is a local functional in time of the intermediate scattering func- 
tions, linear in the tagged-particle intermediate scattering function. 
The coupling coefficients referred to as vertices are determined by 
structural quantities only.^l In addition to the mode-coupling ap- 
proximation we rely on a generalized hydrodynamics approxima- 
tion, i.e. we replace the memory kernel by its long-wavelength 
counterpart M(q,t) n> M(0,t) =: m(t). The microscopic expres- 
sion is given explicitly by 



m(t) 



(5) 



For independently distributed scatterers the direct correlation func- 
tion between matrix and tracer is provided by the spatial Fourier 
transform of the Mayer function.^ For hard exclusion one obtains 
c s (k) = —47ccr 2 }i(kcr)/k, where jj(-) denotes a spherical Bessel 
function. The critical behavior encoded in this theory coincides 
with the original self-consistent kinetic theory. mmiThe mathemat- 
ical structure of the simplified mode-coupling theory is identical to 



the one of the mode-coupling theory of mixtures for which various 
properties have been proven rigorously, example, for Brown- 
ian dynamics there exist unique solutions with positive spectra, i.e. 
correlation functions, which are completely monotone functions. 

The MCT approach for the Lorentz model is encoded also in 
Rrakoviack's MCT description^ for partially pinned fluids in the 
limit of a dilute fluid component. There the interaction with the 
disordered matrix is incorporated on a more sophisticated level 
than in our equations and enters via connected and disconnected 
(blocked) parts of the direct correlation function. Our equations 
coincide with his approach within the generalized hydrodynamics 
approximation, 8 provided the disconnected correlation function is 
ignored, which is a common approximation in integral equation 
theories. 

The predictions of the theory for the localization transition can 
be discussed most conveniently in the Fourier-Laplace domain, 
convention F(q,z) = i Jo'F(q,t)e lzt dt. Then Zwanzig's equation 
of motion yields 



(6) 



Z + iv s + m(z) 

By linearity m(z) is obtained by the same formula as eq. |5]( with 
the intermediate scattering functions replaced by their Fourier- 
Laplace transform F(q,z). Combining both equations leads to sin- 
gle non-linear equation for the unknown kernel m(z). Abbreviating 

fi(z) = z + ivj + m(z), one has to solve for 

iZav . r z 2 /i( z ) 2 

a{k) 2 -zfi{z) 

(7) 

A bifurcation occurs at a critical density n^ CT = 9/47rcr 3 where 
1 = {nf CT /3) J A7ik 2 &k(2n)- 3 c s (k) 2 . Directly at this critical den- 
sity zfi(z) becomes small for frequencies z — > 0, and upon expand- 
ing one obtains 



-/ i 2 ■ n f°° 47rk 2 dk .,,7 f , 
zn{z)-z -iv s z=- j o -j2^T c ^ k ) + 



-z 2 -iv s z = 



„MCT r4nk 2 dk 



< Xk) 2Z ^+ff(zm 512 (8) 



3 Jo (2tt) 3 n ' Q(*)2 

where the correction term arises from the singular behavior of 
the denominator for long wavelength. The integral can be per- 
formed yielding — izv s = (6cr 2 /5v 2 )[— zm(z)} 2 as leading contri- 
bution to the equation. Thus at the critical point the memory ker- 
nel displays singular behavior —zih(z) = (— izf. s ) 1 ^ 2 v 2 /3cr 2 + &{z) 
with t s = 135<t/8v as has been calculated already by Gbtze et 
a /J 31 | 32 | xhis implies for the intermediate scattering function for 
long-wavelengths and low frequencies 



-zF(q,z) 



l+q 2 0- 2 (-izt s )- l l 2 



(9) 



The right-hand side displays spatio-temporal scaling behavior, it 
depends only on the combination zlq" which implies a walk di- 
mension rfJ^ CT = 4. Furthermore for large wavenumbers qcr S> 
(|z|f,) _ly ' 4 , the scaling function decays as q~ 2 from which we in- 
fer the fractal dimension d^ CT = 2. Transforming back to the tem- 
poral domain, we infer that at the critical density, F(q,t) displays 
scaling behavior for long times and small wavenumbers similar to 
eq. |3| with exponents that are simple numbers. 
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The precise value of the critical density within mode-coupling 
theory is determined by the properties of the vertex, i.e. the cou- 
pling coeffients. In the standard factorization procedure used here, 
«J' ICT cr 5 = 9/4n ~ 0.716 is surprisingly close to the percolation 
threshold n c cr 3 ~ 0.838. The situation is less favorable in two di- 
mensions, where MCT yields a value of 2/n ~ 0.637 whereas the 
percolation occurs at 0.359. Interestingly, if one extrapolates the 
rigorous low-density expansion for the case of Brownian tracers,^ 
the situation is reversed. There the low-density expansion leads 
to 3/2n ~ 0.48 in three dimension whereas in two dimensions one 
arrives at l/n « 0.3 1 8. 

Strictly speaking the scaling properties discussed above are a 
consequence of the generalized hydrodynamics approximation. 
Using the full wavevector-dependent theory these exponents 
hold only in an intermediate regime. This phenomenon of spuri- 
ous long-wavelength singularities in MCT approaches is now well- 
understoodEDand we do not repeat the discussion here. Ultimately 
it remains unclear how to properly handle the emergence of frac- 
tal structures and the divergent length scale^nevertheless we shall 
show that certain aspects predicted by the theory have a close coun- 
terpart in the simulation results. 

Rather than relying on asymptotic analysis we solved the mode- 
coupling equations numerically in the temporal domain. The in- 
tegral in the memory kernel is one-dimensional in the wavenum- 
ber k when expressed in spherical coordinates. To account for the 
long-wavelength behavior we have used a wavenumber discretiza- 
tion with N s logarithmically-spaced wavenumbers from a smallest 
fc miti up to a wavenumber Ak. From there we use an equidistant 
grid with spacing Ak up to a largest wavenumber fc max . We used 
N s = 50, k min = 10~*/R, Ak = 0A/R and fe max = 24/R. With these 
parameters, the critical density is slightly shifted to n* ~ 0.748 but 
converges to the exact value on finer discretizations. The singular 
behavior of the memory kernel given above matches the numerical 
results to within 1% of its amplitude at long times. 

A standard MCT algorithm was used to solve the equation on an 
equidistant time grid, which is repeatedly coarsened by factors of 

2 in order to be able to calculate the solution over many decades 
in time. The memory kernel was numerically integrated on the 
wavenumber grid with the trapezoidal rule. 

We have checked that the spurious long-wavelength singularities 
do not introduce a qualitative difference in the ranges displayed in 
the figures, with the exception of the mean-square displacement 
(not shown), by comparing to the full wavenumber-dependent the- 
ory. For a detailed discussion of the implications of using a dis- 
cretized wavenumber-dependent mode-coupling functional we re- 
fer to Schnyder et alJ^SI 

3 Simulation 

We generate trajectories by event-driven molecular dynamics sim- 
ulations at the critical density for system sizes L = 200<x corre- 
sponding to 6,704,000 obstacles. We impose periodic boundary 
conditions to minimize finite-size effects. Initial particle posi- 
tions are sampled only on the remaining infinite cluster, which 
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Fig. 1: Intermediate scattering function F ao (q,t) for tracers exploring 
only the percolating cluster. Time is measured in units of t = tr/v. The 
simulation is performed directly at the critical density n* = 0.838 and com- 
pared to the mode coupling theory. The dotted black line corresponds to 
an exponential relaxation for comparison. The full black line is a fit to a 
stretched exponential exp(— (t/r q )P), fi = 0.39. 

we identify by a numerical Voronoi tessellation]^ Out of 96 re- 
alizations of the disorder only 44 displayed an infinite cluster, con- 
sistent with the theoretical expectation that in finite systems the 
transition is smeared. We extract a critical density 49 n* = 0.838 
which is slightly lower than the value 0.839 found for L = 200cr in 
Ref s .|26|27|The difference can be rationalized by the bias to obtain 
half of the samples to be percolating. 

Correlation functions are obtained by moving-time averages for 
lag times up to 3 • 10 9 r o using a standard blocking scheme 50 where 
the trajectory is sampled essentially on a logarithmic time grid. 
Our results are extracted from 7 independent trajectories in each 
of the 44 percolating configurations. We measure the incoherent 
intermediate scattering function directly as defined in eq. |2j with 
the constraint that the wave vectors are consistent with the periodic 
boundaries. The smallest wavenumber that can be achieved in the 
simulation is limited by the system size to 2n/L, i.e. q = 0.03 l<x~ 1 . 

The main computational effort lies in the identification of the 
percolating cluster and to locate the critical density with sufficient 
accuracy. The limitation of the system size is set by computer 
memory necessary (25 GB) to analyze the Voronoi network. Once 
the infinite cluster is identified, a single trajectory requires typi- 
cally 12 hours on a prevalent 3 GHz CPU, resulting in 150 days of 
computing time per system size and density. 

4 Results and discussion 

Simulation results for the three-dimensional Lorentz model have 
been presented mostly in terms of the mean-square displace- 
ment |26|27|30] Spatio-temporal information beyond the second mo- 
ment is encoded in the intermediate scattering function, as shown 
in fig. [T] for the obstacle density directly at the transition point. 
The time window in the simulation covers 9 non-trivial decades 
and the data displays practically no noise. The range of wavenum- 
bers includes a factor of almost 200 corresponding to the smallest 



*We have employed the free voro++ package originally developed by Chris 
Rycroft, see http : / /math . lbl . gov/voro + + /. 
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Fig. 2: Relaxation time T q of the simulated intermediate scattering func- 
tion at the critical density n* = 0.838 for different wave vectors q and sys- 
tem sizes L. The thick grey lines indicate a power law expected from the 
scaling hypothesis and serves as a guide to the eye. Inset: Rectification plot 
of the relaxation times testing the asymptotic power-law behavior. 



Fig. 3: Intermediate scattering function at the critical density, double- 
logarithmic representation. The terminal relaxation approaches a power 
law for small wavenumbers. As guide to the eye: solid black lines are 
power laws r^r/^ ana ; f~'/ 2 respectively. 



length scale of an obstacle diameter cr to the box size L = 200cr. 
A pronounced bump becomes apparent at the collision time t c ~ t 
for wavenumbers probing the microscopic scales, familiar from the 
dynamics of dense liquids For long wavelengths the decay be- 
comes considerably non-exponential, rather the shape is charac- 
terized by significant stretching. A fit to a Kohlrausch- Williams- 
Watts function exp(— (t/Tqf) provides a nice description of the 
data down to a value of 0.4 and yields a rather low stretching ex- 
ponent /3 = 0.39. For even lower values of q, the stretching be- 
comes even more significant. Within the long wavelength regime 
the shape of the ISF become insensitive to changes in q, only the 
characteristic time scale varies. The MCT results for the critical 
obstacle density n^ CT are included in the figure for comparison. 
For moderately small wavenumbers, 0.4 < q < 2.5, the agreement 
is reasonable, although a pronounced tail appears in the simulation 
data that is not captured by the theory. For lower wavenumbers the 
shape of the relaxation curves is still similar to the simulation data, 
yet the time scales start to deviate significantly. 

Decay times r q defined by Foo(q,Tq) = 1/e vary over 8 orders 
of magnitude. Note that for ordinary diffusion in the hydrody- 
namic regime, one would expect a change of r q by only 4 or- 
ders of magnitude, since then F(q,t) = exp(—Dq 2 t) and T q oc q~ 2 . 
The relaxation times r q as a function of q are exhibited in fig. [2] 
using a double-logarithmic representation. The data approach a 
power law T q ~ q~ d ™ in the regime of small wavenumbers and, 
in principle, allow the determination of the walk dimension d w . 
The rectification plot (inset of fig.[2]( illustrates that the power-law 
behavior is approached only asymptotically. We fix the scale fac- 
tor in T q /t ~ A T {qa)^ dvi to A T = 65 ± 10. From the slope one 
can in principle obtain the value for d w which is compatible with 
our earlier results 2930 extracted from mean-square displacement 
SrL(t) ~ f 2/,rfw . Of course, d w can be determined more accurately 
from drt(t) since the time window spans 8 decades, whereas the ar- 
range is varied only by a factor of 200, hence we rely on our earlier 
estimate in the remainder of this article. 

To investigate the behavior at long times we display the sim- 
ulation data in a double-logarithmic representation in fig. [3] In 



the long-wavelength regime the curves approach a power-law re- 
laxation for t — > oo. From the scaling theory at the critical point, 
eq. (BJ, for small rescaled arguments one anticipates Fao(q,t) ~ 
t -d,JTj _j, ^ wnere ^ s _ 2df/d w is known as the spectral dimen- 
sion.^ Our data are compatible with the value d s = 1.054 ob- 
tained from the previously determined numbers for df and d w . The 
spectral dimension for three-dimensional Lorentz models is sur- 
prisingly close to the MCT result d^ 1 = 1, in fact based on the 
window accessible in today's computer simulations these values 
are indistinguishable. 

In addition to reading off the critical exponents from relaxation 
times T q and the power-law terminal relaxation, we investigate the 
scaling properties, see fig.|4] Shifting the data according to rescaled 
time t ~ q d "t the curves almost superimpose. Upon decreasing the 
wave number, the data collapse improves and asymptotically a uni- 
versal master curve is approached. Deviations at moderately small 
wave numbers are predominantly due to the poor convergence of 
Tqlq~^" (fig. [2J whereas the shape appears to be approached more 
rapidly. In fact using the measured T q rather than its asymptotic 
prediction enforces data collapse at the point ( 1 , e~ 1 ) and the scal- 
ing function is nicely followed. The scaling plot is sensitive to fi- 
nite size corrections and fine-tuning to the critical density, as shown 
in the top panels of fig. [4] Interestingly, MCT displays an approx- 
imate data collapse for small wavenumbers, as is explained by the 
small difference in d w . 

The MCT estimate for the walk dimension is d™ CT = 4 is suffi- 
ciently close to the measured value d w =4.81 to analyze the data 
with a rescaled time f ~ q^t. Such a rescaling is displayed in 
fig. [5] and indeed the simulation data show reasonable data col- 
lapse, given that the original data exhibit a spread of six decades. 
As expected, the numerical solutions of the MCT equations now 
collapse perfectly for small q. 

Since the scaling prediction is an asymptotic law, valid for long 
times and small wave numbers, an even more sensitive test is ob- 
tained by plotting t d '^ dv ' F„(q,t) vs. rescaled time ?~ q dw t. Then 
the long-time behavior of F„ is highlighted and the approach to- 
wards a power-law relaxation is mapped to a finite limit. Our simu- 
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ltr 2 10" 

i= A- 1 (qG) d 't/t 

Fig. 4: ISF rescaled as a function of rescaled time f = A™ 1 (gcr)' 1 '" t / t a for 
small wavenumbers <j. The lower panel corresponds to the largest system 
size investigated (L = 200cr) at the critical density (n* = 0.838). The top 
panels are for smaller system sizes (left, middle), or larger distance ((n c * 
—n*)ln* = 0.01) to the critical density. 
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Fig. 5: For comparison: rescaling of the intermediate scattering function 
with rescaled time (q<r) 4 t/t„ as suggested by MCT for the same data as in 
Fig.[l] 



lation data (see fig.[6]( nicely follow the expected scaling for small 
f; for long times, however, deviations become visible. Neverthe- 
less, one infers slow convergence for small wavenumbers q — > 0. 
In principle, the scaling hypothesis can be extended to include uni- 
versal corrections to scalingESIto rationalize the residual spread. 



1.0 



o.s 



^.0.6 
4 0.4 



0.2 



0.0 



MCT results: 

q = 1.20 

q = 0.40 

q = 0.28 

q = 0.14 

q = 0.07 




Simulations: 

5 = 1.13 

g = 0.50 

g = 0.25 

g = 0.13 

g = 0.06 



nr 2 io° 

t = A-^qa^t/to 

Fig. 6: Rescaled and rectification highlights the approach to a master 
curve. The scaling variable for the simulation data is t = A^ l (tlt )(q(rY v ' 
and spectral dimension d s = 1.054, respectively (in the inset) t = 
(r/f )(^o") 4 and spectral dimension 1. 



The success of the theory is encouraging on the one hand, but 
delusive upon closer inspection. First, we chose to compare only 
the motion on the infinite cluster to MCT, rather than the all-cluster 
average dynamics. While a dynamic scaling hypothesis can be 
worked out based on the percolation picture, see e.g. Refs. 5 1 and 
l46l the exponents start to deviate significantly from the measured 
ones. Second, there is a conceptual problem, since MCT does not 
distinguish between the two types of tracers. While the tracers 
on the infinite cluster are ergodic, compatible with the MCT no- 
tion, the all-cluster average motion is not, since some particles will 
always be trapped in finite pockets of the host structure. These 
holes are present at arbitrarily low obstacle density, leading to 
non-vanishing long-time limits of the ISF sensitively dependent of 
wavenumber. 29 This peculiarity supports the idea that we have pur- 
sued in the current paper to make the comparison with the infinite 
cluster motion only. Yet, for obstacle densities above the transition, 
MCT predicts a continuous increase of the long-time limits of the 
ISF. Since the infinite cluster ceases to exist a comparison here is 
excluded. 



5 Summary and Conclusion 

We have measured the intermediate scattering function for a three- 
dimensional Lorentz model directly at the percolation transition for 
particles confined to the infinite cluster. The dynamics becomes 
anomalous due to the underlying spatial fractal and a scaling hy- 
pothesis is expected to hold. We have tested various aspects of 
spatio-temporal transport thereby corroborating the notion of self- 
similarity beyond the second moment. In particular, we find a 
power-law relaxation in the long-wavelength regime strikingly dif- 
ferent from the two-step relaxation scenario known from glasses 
and supercooled liquids. 5 We have compared our simulation re- 
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suits within a simplified mode-coupling theory that reproduces the 
critical properties of a self-consistent mode-coupling kinetic the- 
ory exponents predicted by this approach are sufficiently 
close to the measured values to allow for a quantitative test of the 
theory and surprisingly good agreement is found. 

While the agreement with MCT should not be overstated for 
the reasons discussed above, it does highlight that the theory con- 
tinues to provide qualitative, and sometimes semi-quantitative, 
predictions also in the vicinity of certain percolation-type transi- 
tions. This is of importance, since there is now emerging a re- 
newed interest in anomalous transport and its interplay with glass- 
transition phenomena, for example in the case of systems with self- 
generated disorder (such as size-disparate binary mixtures). Since 
true asymptotic scaling behavior emerges only rather close to the 
transition and often outside experimentally accessible windows, 
MCT proves useful as a theory that can also, beyond scaling anal- 
ysis, predict generic (non-universal, but typical) behavior. 

While we have clarified what the intermediate scattering func- 
tions should look like, it remains a challenge for the future to derive 
a mode-coupling theory from first principles that accounts for the 
percolation transition. A possible route would be to build the the- 
ory on the obstacle clusters that form and use renormalized vertices 
where the geometry of the obstacle parcours is reflected properly. 
The idea for such a cluster-MCT has been promoted for the prob- 
lem of weak gelation of attractive colloids, 52 where colloids stick 
to clusters moving together as entities. However, it remains open 
how such a theory can be achieved in practice. A second promising 
approach is to employ a diagrammatic expansion 53 for the dynam- 
ics to first rederive the simplified MCT equations and then identify 
what processes have to be included to account for the underlying 
fractal geometry. 
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